Shannon dimensionality of quantum channels and its application to photon 

entanglement 
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We introduce the concept of Shannon dimensionahty as a new way to quantify bipartite entan- 
glement as measured in an experiment. This is apphed to orbital- angular- momentum entanglement 
of two photons, using two state analyzers composed of a rotatable angular-sector phase plate that 
is lens-coupled to a single-mode fiber. We can deduce the value of D directly from the observed 
two-photon coincidence fringe. In our experiment, D varies between 2 and 6, depending on the 
experimental conditions. We predict how the Shannon dimensionality evolves when the number of 
angular sectors imprinted in the phase plate is increased and anticipate that D :^ 50 is experimen- 
tally within reach. 
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Photons can be entangled in various degrees of free- 
dom. The most extensively studied variety involves the 
polarization degrees of freedom, of which there are in- 
herently two per photon. In a typical EPR-Bell type 
experiment, the state analyzers are polarizers, and when 
their relative orientation is scanned, this gives rise to a 
sinusoidal coincidence fringe [1]. This particular shape is 
characteristic of the two-dimensional nature of polariza- 
tion entanglement. 

Recently, much attention has been drawn to bi- 
partite entanglement involving more than two degrees 
of freedom. With increasing dimensionality, quantum 
entanglement becomes correspondingly richer. High- 
dimensional entanglement is predicted to violate locality 
more strongly and to show more resilience to noise P, W\. 
From an applications perspective, it holds promise for 
implementing larger alphabets in quantum information, 
e.g. quantum cryptography [^L and for an increased se- 
curity against eavesdropping Q. High-dimensional en- 
tanglement can be studied employing the frequency-time 
[g.] or position-momentum degrees of freedom, the latter 
fmving been demonstrated for both the transverse linear 
and orbital-angular-momentum degrees of freedom 

It is crucial to have a quantifier of the dimensional- 
ity of entanglement as measured in an experiment [11]. 
In this Letter we introduce such a quantifier, using con- 
cepts from classical information theory in the spirit of 
Shannon [12^. We apply these ideas to orbital-angular- 
momentum entanglement, inserting appropriate angular 
state analyzers in the beamlines of a parametric down- 
conversion setup. We have realized a Shannon dimen- 
sionality 2 < D < 6 and we argue that D :^ 50 is within 
reach. 

In classical information theory [12], the number of in- 
dependent communication channels of a signal is known 




as the Shannon number. The signal being the state of 
a physical system, the Shannon number is also referred 
to as the number of degrees of freedom, or the number 
of modes, of that system [13,, ^4]. For example, a signal 
encoded in the polarization degrees of freedom of a light 
beam has a Shannon number equal to 2. 

When dealing with a bipartite quantum system in an 
entangled pure state |?/^) G J(f = J(fA ^ the usual 
measure of the effective dimensionality of the Hilbert 
space in which the state lives, is given by the Schmidt 
number K |l5| 



K 



Tr^(p^) Txb{p\)' 



(1) 



Here, pA = TrB(|V^)(V^|) and pB = TrA(|V^)(V^|), are the 
reduced density matrices representing the states of the 
two sub-systems A G Ji^A and B G Ji^B^ respectively. 
Although a system may have infinitely many degrees of 
freedom, any actual measurement apparatus has effective 
access only to a finite number of them, say D. Such a 
dimensionality D is referred to as the Shannon number 
of the measurement apparatus. 

Consider an experiment measuring correlations be- 
tween the two subsystems A and B. There are two 
measuring apparatuses, say ^a(<^) and ^b{P)^ inter- 
acting with subsystems A and 5, respectively, where a 
and P label possible settings of the two apparatuses. For 
a given setting ^ G {a, detector represented 
by the projection operator r(^) = where 
X G {A, 5}, and \X{^)) is the state in which the system 
X is left after measurement. 

If a Von Neumann- type projective measurement is per- 
formed, the set of states {|X(^))}^ obtained by varying 
^ is complete and orthonormal, namely 



{x{o\x{e)) = s, 



1, 



(2) 
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where the measurement operators f (^) are Hermitean 
and idempotent. The number of these operators is equal 
to the dimension of the Hilbert space of the measured 
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quantum system [l6]. However, in many situations non- 
orthogonal measurements are made and Eqs. ([2]) do not 
hold [13 • In this case, the number of projection operators 
r(^) does not give the dimension of the Hilbert space 
of the measured system, and a new criterion must be 
introduced. 

Let us therefore consider finite-dimensional systems, 
say dim(.;^) = L, and rewrite Eq. ([2j) for the case of 
non-orthogonal measurements as 

(x(0|x(O)=^ce, Er(0=7, (3) 

where G = [g^^'] is a matrix of size L x and 7 is an 
Hermitean operator. The eigenvalues of 7 give the 
detector's 'sensitivity' to the corresponding eigenmodes. 
In general, a detector will not be equally sensitive to all 
eigenmodes and some are substantially larger than oth- 
ers. The effective dimensionality D < L of the Hilbert 
space ^ where the measured system lives can be quan- 
tified as the Hilbert- Schmidt norm of the eigenvalue dis- 
tribution [m 

This dimensionality should be interpreted as the effective 
Shannon number of information channels [12, 

The isomorphism of Eq. ([T]) and Eq. (|4]) suggests a re- 
lation between the Schmidt number K and the Shannon 
dimensionality D. The nature of such relation becomes 
clear if one notes that since the operators r(^) are Her- 
mitian and positive semidefinite, the operator 7 may be 
interpreted as a density matrix acting in J^x 0- Thus, 
if we think of 7 as a reduced density matrix of a bipartite 
system, then K and D are formally the same. However, 
it is important to note that while K furnishes the dimen- 
sionality of the generated entanglement, D gives the ef- 
fective dimensionality of the space ^ that can potentially 
be probed and it is a property of the projection apparatus 
only. The dimensionality of the measured entanglement 
is a joint property of the generated system and analyzers, 
but simply amounts to D as long as ^ D ^. 

Next, we apply our formal theory to an experiment 
on orbital-angular-momentum entanglement of two pho- 
tons, in order to illustrate how detector characteristics 
bound the measured entanglement to an effective Shan- 
non dimensionality while probing a generated state 
with Schmidt number K ^ D (and ^ D ^). Our exper- 
imental setup is depicted in Figure [TJ Pumping a BBO 
non-linear crystal with a 150 mW Kr+ laser beam at 
A = 413 nm, we produce spatially entangled photons by 
means of spontaneous parametric down conversion. The 
state we generate is of the form |^) = X1^V%^I0I ~ 0^ 
where |/) denotes the orbital angular momentum eigen- 
mode of order /: {(j)\l) = exp(i/0)/v/27r, with (/) the az- 
imuthal angle [2^. Employing Type-I collinear phase 
matching, we collect the full emission cone and with the 
experimental parameters of our setup (beam half-waist 



at the position of the crystal Wq = 250 /im and crystal 
length 1 mm) we obtain an azimuthal Schmidt number 
K 31 [m, [13, [23. The twin photons are spatially 
separated by means of a non-polarizing beam splitter. 

Each arm of the setup contains an angular state ana- 
lyzer, composed of an angular phase plate that is lens- 
coupled to a single- mode fiber (see Fig. [1]) [24]. The 
angular phase plates carry a purely azimuthal variation 
of the optical thickness. As in polarization entanglement 
[H , the phase plates are rotated around their normals and 
the photon coincidence rate is recorded as a function of 
their independent orientations |lOl |. 

The combined detection state of the two angular- 
phase-plate analyzers, each acting locally, can be ex- 
pressed as 

{5} 

where a and P denote the orientations of the two phase 
plates, respectively [25]. The complex expansion coef- 
ficients are fixed by the physical profile of the an- 
gular phase plate and obey the normalization condition 
1 7/ 1 = 1. In general, the detection state constitutes a 
non-uniform superposition of orbital angular momentum 
modes. When the angular phase plates are rotated over 
a or all modes in the superposition rephase with re- 
spect to each other, yielding a set of detection states of 
the type Eq. (|3]). The effective Shannon dimensionality 
that is so probed is given by Eq. (jlj). It is the average 
number of modes captured by an analyzer when its phase 
plate is rotated over 360°. 

As we have recently shown in Ref. [18] , the Shannon di- 
mensionality is straightforwardly deduced from the shape 
of the experimental coincidence curve; it is the inverse 
of the area underneath the peak-normalized coincidence 




FIG. 1: (color online). Experimental setup. Orbital-angular- 
momentum entangled photons are emitted at 826 nm by a 
BBO crystal, cut for Type-I collinear phase matching. A 
thin GaP wafer serves to eliminate the pump beam. The 
two-photon field can be clipped with an aperture. The twin- 
photons are spatially separated by a beam splitter and imaged 
on the angular phase plates (/2 = 4/i = 40 cm). Just be- 
hind the phase plates, the frequency-degenerate photons are 
selected by interference filters (not shown), centered around 
826 nm with a 10 nm width. The phase plates (shown are 
quarter- sector plates) are oriented at angles a and /3, and 
photon counts are rendered by a coincidence circuit. 
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fringe, obtained when rotating one of the phase plates. 

In our experiment, we have used angular-sector phase 
plates; these have a single arc sector, characterized by 
the angle ^, whose optical thickness is A/2 greater than 
that of the remainder of the plate [24]. The part of the 
field that crosses this sector thus flips sign. The phase 
plates are manufactured from fused-quartz plane-parallel 
plates, having a wedge angle of 0.25" . They are processed 
by a combination of photolithography, wet etching, depo- 
sition and lift-off, resulting in a well-defined mesa struc- 
ture, with a transition region that is typically 20 jam 
wide. The insets of Figure [2] show two such plates; a 
half-sector plate {6 = tt) consisting of two equal halves 
that are phase shifted by tt; and a quarter-sector plate 
{S = 7t/2) having one quadrant 7r-phase shifted with re- 
spect to the remainder of the plate. 

For state analyzers that are equipped with such sector 
phase plates, the Shannon dimensionality is given by [Ts*] 



I)(27r-(5), 



(I) 



, Se[0,7r], 
5 e [7r,27r]. 

(6) 

For (5 = we find the trivial result D = 1; a planar 
plate does nothing. For ^ = tt, i.e. a state analyzer 
equipped with a half-sector plate, we arrive at = 3. 
For an analyzer equipped with a quarter-sector plate we 
find D = 6. This is the maximum value for a single 
angular-sector phase plate. We note that for our setup 
indeed K:^ D. 

In the experiment, we scan one angular-sector phase 
plate over a 360° rotation, the other remaining fixed, 
and measure the coincidence rate. In terms of Klyshko's 
picture of advanced waves [2^, valid when K ^ the 
resulting shape of the coincidence curves can be explained 
in terms of the mode overlap of the two state analyzers. 
Figure[2](a) shows experimental results obtained with two 
half-sector plates = tt), having step height of 0.48 A. 
The data points form a double parabolic fringe, consis- 
tent with theory (solid curve). The maxima at 0° and 
180° are sharply peaked. The zeros of the fringe are very 
deep; less than 10 counts per 10 seconds. The maximum 
coincidence rate is of the order of 6.5 x 10^ per 10 sec- 
onds, compared to 10^ single counts. We verified that 
the coincidence rate depends on the relative orientation 
between the two phase plates only, the fringe visibility 
being > 99% for all cases studied. This basis indepen- 
dence is the key aspect of quantum entanglement. From 
the area underneath the data we deduce the experimen- 
tal value D = 3.0. Note that a parabolic fringe was also 
reported in Ref. [To^, obtained with non-integer spiral 
phase plates. We conclude that also in that case D = 3. 

An aperture, positioned inside the telescope, allows us 
to control the number of detected modes (see Figure [T]). 
Because of the anti-symmetric profile of the half-sector 
plate, the detection state contains only odd expansion 
terms (see Eq. (|5])) in a fashion ji = When the 

aperture size is reduced, higher-order orbital-angular- 
momentum modes are cut off so that, eventually, only 



the modes / = 1 and / = — 1 survive. We then expect 
a sinusoidal fringe, analogous to two-dimensional polar- 
ization entanglement [l|. In the experiment, we observe 
that the coincidence curve is gradually transformed from 
parabolic to sinusoidal when the aperture gets smaller. 
Using an aperture of 600 /im diameter, we are in an in- 
termediate regime (squares, D = 2.1), while using a 400 
/im diaphragm yields a curve that resembles a sine very 
well (triangles, D = 2.0). The dashed and dotted curve 
are theoretical predictions. 

To achieve I) = 6, we use two quarter-sector plates 
{S = 7r/2), carrying an edge discontinuity deviating less 
than 3% from A/2. The circles in Figure [2fb) show our 
experimental results, revealing a coincidence curve which 
is parabolic for |a — /3| < 90° and equal to zero otherwise, 
in agreement with theory (solid curve). We find D = 
5.8, in very good agreement with the expected value of 6 
mentioned above. 

The maximum value of the Shannon dimensionality 
that can be achieved with a phase plate having but a 
single sector is D = 6. Can one reach higher values of D 
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FIG. 2: (color online). Coincidence count rate vs. the relative 
orientation of the two state analyzers. Points denote experi- 
mental data, the curves are theoretical predictions, (a) Half- 
sector plate. The parabolic fringe (circles) is a signature of a 
dimensionality larger than two: we find D — 3.0. Truncat- 
ing the number of modes, by closing the aperture, gradually 
reduces the parabola into a sine of dimensionality 2.0 (tri- 
angles), (b) Quarter- sector plate. The piece- wise parabolic 
fringe yields an experimental dimensionality of 5.8 (circles), 
where theory predicts D = 6. 
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by using plates with more sectors? To answer this ques- 
tion, we consider plates with N sectors that are phase 
shifted by tt with respect to interjacent regions. For each 
choice of sector angles, we calculate the expansion coeffi- 
cients and, subsequently, D (see Eq. (j4]) and ([5])). 
Next, we maximize D by adjusting the sector angles us- 
ing a Monte-Carlo random-search algorithm. The result 
is plotted in Figure [3l showing a graph of the maximum 
value of D versus the number of mesas N. For 10 such 
sectors, we find D = 49.9. The insets show the optimal 
phase plates for N = 1 (quarter-sector plate). A/" = 2, 
and = 3. 

In conclusion, we have introduced the effective Shan- 
non dimensionality as a novel quantifier of entanglement 
as measured in an actual experiment. We have demon- 
strated its significance to the case of two-photon orbital- 
angular- momentum entanglement. Using angular-sector 
phase plates, we have achieved Shannon dimensionalities 
up to I) = 6. We anticipate that it is feasible to probe 
dimensionalities as high as 50, using multi-sector phase 
plates. These can be manufactured by means of photo- 
or e-beam lithography as in diffractive-optics technology. 
Alternatively, the use of adaptive optical devices, such 
as spatial light modulators or micro- mirror arrays, seems 
promising, particularly because of their versatility with 
regard to plate patterns. However, the ultimate limit to 



the Shannon dimensionality is constrained by the angular 
Schmidt number of the source; using periodically poled 
crystals, such as PPKTP, K ~ 100 is viable for realistic 
values of pump-beam waist and crystal length, without 
loss of count rates [27]. 

This work was supported by the Stichting voor Funda- 
menteel Onderzoek der Materie. 




FIG. 3: (color online). Maximum dimensionality that can 
be accessed with sector phase plates having 2N angular sec- 
tors alternatingly phase shifted by tt. The insets show the 
optimized plates for N = 1, N = 2 and N = 3. 
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